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Abstract 

A QCD analysis of the world data on polarized deep inelastic scattering is presented in 
leading and next-to-leading order. New parameterizations are derived for the quark and 
gluon distributions for the kinematic range x e [10~ 9 ,1], Q 2 e [1, 10 6 ] GeV 2 . The ex- 
trapolation far outside the domain of the current measurements is given both to allow 
for applications at higher values of Q 2 and to be able to calculate integral properties of 
the present distributions. The values of Aqcd and a s (M z ) are determined. Emphasis is 
put on the derivation of the fully correlated la error bands for these distributions, which 
are also given in terms of parameterizations and are directly applicable to determine ex- 
perimental errors of other polarized observables. The impact of the variation of both the 
renormalization and factorization scales on the value of a s is studied. Finally we perform 
a factorization-scheme invariant QCD analysis based on the observables g\{x,Q 2 ) and 
dgi(x,Q 2 )/dlog(Q 2 ) in next-to-leading order, which is compared to the standard analysis. 
A series of low moments of parton densities, accounting for error correlation, are given to 
allow for comparison with results from lattice simulations. 



1 Introduction 



The nature of the short-distance structure of polarized nucleons is one of the central questions of 
present day hadron physics. Nucleons are composite fermions and their spin should be obtained 
as a superposition of the spins and angular momenta of its constituents, the quarks and gluons. 
At large enough space-like 4-momentum transfers — q 2 = Q 2 and large energy transfer v = 
2p.q/M, with p and M the nucleon 4-momentum and mass, the QCD-improved parton model 
applies. A QCD analysis at leading twist level allows to determine the polarized parton densities 
and the QCD scale A QCD up to next-to-leading order (NLO) [1-3]. The distribution functions 
derived can be used to calculate other leading twist hard scattering cross sections of polarized 
nucleons. 

During the recent years the remarkable growth of deep inelastic scattering data off polarized 
targets [4-13] allows to perform a detailed analysis. Here the different systematic effects of 
the data have to be taken into account. Previous analyzes [14-16] [] present parameterizations 
only for the central values of the parton distributions. It is, however, desirable to derive at 
the same time parameterizations of the polarized parton densities and those of their la error, 
taking into account all error correlations being available. With these parameterizations at hand 
it becomes possible to give estimates on the uncertainty of the measurement of other hard 
scattering processes of polarized targets w.r.t. the knowledge of the parton densities. In the 
analysis the polarized parton distributions are first determined at a reference scale Q$. The 
aforementioned parameterizations are found evolving both the central values and their errors 
from this scale to all scales Q 2 of interest, which we choose as 1 GeV 2 < Q 2 < 10 6 GeV 2 for 
x = Q 2 /2p.q e [10~ 9 , 1], and are given in form of fast and accurate grid interpolations. 

Unlike the case in most of the previous analyzes the QCD analysis performed in the present 
paper determines the Q CD-scale Aqcd along with the parameters of the parton densities at 
the initial scale Qq. Moreover, positivity constraints are not imposed from the beginning but 
first left to the fit. A systematic investigation of the x 2 ~P r °fil es during the analysis showed, 
which of the parameters of the parton densities have enough sensitivity to be measured from 
the current polarized deep-inelastic scattering world data. In this context it turns out that 
the current inclusive data do not yet allow to pin down the flavor structure of the sea quarks. 
Therefore we are going to use flavor SU(3) in the current analysis and leave a refinement to 
a further analysis including as well semi-inclusive data.f] As the subsequent analysis shows, 
the small-x behavior of the parton densities cannot be reliably fixed using the current data. 
Special dynamical effects [[|], |2]J may be present at smaller values of x < 5 • 10 -3 than currently 
probed, however, reliable estimates cannot be given at present due to the lack of calculations 
on higher order resummed corrections in this kinematic range. Therefore we present two sets of 
parameterizations which could be both accommodated to the present data. 

The standard analysis is performed using the MS factorization and renormalization scheme. 
In this way two independent scales, fi 2 and fi 2 R , are introduced into the parton densities. The 
observables, g^x.Q 2 ) and g™(x,Q 2 ), are independent of these scales which leads to two asso- 
ciated (matrix-valued) renormalization group equations for the Wilson-coefficients and parton 
densities, which are the evolution equations. Referring to also other observables, as e.g. the 
slopes dgi' n (x,Q 2 )/d\ogQ 2 , it is possible to each order in the coupling constant a s = a s /(47r) 
to formulate physical evolution equations |2], |23| of observables itself, which removes the depen- 
dence on fi 2 . This method of analysis is conceptually of high interest in future analyzes, having 

1 In |l7| another factorization scheme is referred to. For a list of older parameterizations see [ fl8| . 
2 The latter analysis can only be performed if the relevant fragmentation functions are known, cf. [pi. 
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precise measurements of the respective observables at the scale Qq at hand. Then the only pa- 
rameter to be fitted is Aqcd- We also perform this factorization-scheme independent analysis 
in next-to- leading order to check the stability of the fitted values of Aq CD . 

Currently the measurement of also the moments of polarized parton densities on the lattice 
become more precise [23-25] . Yet a series of systematic and algorithmic effects has to be studied 
in detail, but one may expect that direct comparison between the moments of parton densities 
being extracted out of the world data and the lattice results will be possible in the future. 
This comparison has to include the proper treatment of errors both on the side of the lattice 
calculations and the perturbative analysis. This usually includes correlated propagation of the 
experimental errors through the evolution equations. To allow for comparisons with results from 
lattice calculations we provide a set of lowest moments of parton densities and their combinations 
including their errors in the present analysis. 

The paper is organized as follows. In sections 2 and 3 the standard and the factorization- 
scheme independent analyzes to next-to-leading order are described. A summary on the exper- 
imental data used in the present analysis is given in section 4. Section 5 contains the details 
on the parameterizations of the parton densities and that of their la errors. In Section 6 the 
results of the standard and factorization-scheme independent analysis on the parton densities 
and their errors at the scale Qq as well as their evolution throughout the kinematic range of 
the parameterization 1 < Q 2 < 10 6 GeV 2 are presented and the respective values of Aq CD and 
a s (M|) are given and are discussed w.r.t the choice of the renormalization and factorization 
scales. The covariance matrices of the parameters determined in the QCD fits are provided. In 
section 7 a series of values of the lower moments of the parton densities and their errors, suit- 
able for comparisons with results of lattice simulations, are calculated. Section 8 contains the 
conclusions. An appendix provides informations on the numerical parameterization of different 
sets of polarized parton densities and their errors. 



2 Standard Analysis 



The twist-2 contributions to the structure function gi(x,Q 2 ) can be represented in terms of a 
Mellin convolution of the polarized parton densities Afi and the coefficient functions AC- by 



9i(x,Q 2 ) 



1 



3=1 



1 dz 
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(1) 



where Cj denotes the charge of the jth quark flavor and Nf is the number of flavors. The scale \if 
denotes the factorization scale which is introduced to remove the collinear singularities from the 
partonic structure functions. In addition the above quantities are dependent on the renormal- 
ization scale fi r of the strong coupling constant a s (fi 2 ) = g 2 (n 2 )/{lQir 2 ). The structure function 
gi(x, Q 2 ), as a physical observable, is independent of the choice of both scales. The parton 
densities and the coefficient functions are dependent on these scales and obey corresponding 
renormalization group equations. 

The singlet and non-singlet parton densities which occur in Eq. (|J) are expressed by the 
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individual flavor contributions as 

AS (z,$) 



N f 

E 

i=l 



Af qi (z,$)+Af q . (z,fx 2 f 



and 



A qi NS (z,fi}) = Af qi (z,ii 2 f )+Af § 



z,/i 2 



(2) 



(3) 



respectively, where f g . denotes the polarized quark distribution of the zth flavor. 
The running coupling constant is obtained as the solution of 



da s (Hr 



dlog(>2) 

where the coefficients of the /3-function are given by 



(4) 
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(5) 



in the MS-scheme. Here Ca — 3, T F — 1/2, Cp = 4/3 and Nf denotes the number of active 
flavors. In the QCD fit to the data we extract Aq^d and choose Nf = 4 in Eq. (|5]) whereas only 
the three light flavors are used to represent gi(x, Q 2 ), Eq. ([j]). The expression for Aqq D is given 
by 



a qcd = /vexp I -- 



ft. 



' 1 , Pi N 



(6) 



Subsequently we will also compare the values of a s at the mass scale of the Z-boson, a s (M|). This 
is obtained matching the values of a s at the charm- and bottom-quark threshold M c =1.4 GeV, 
M b = 4.5 GeV using Eq. (|) for the value of A^g D . 

The change of the parton densities with respect to the factorization scale \A = Q 2 is described 
by (matrix-valued) renormalization group equations, the evolution equations, which read 



gAgs (x,Q 2 
<91ogQ 2 
d ( A£(x,Q 2 ) 
dlogQ 2 ^ AG(x,Q 2 ) 



Pj s (x,a s )®Af" s {x,Q 



(7) 
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and <S> the Mellin convolution 

[A <g> B](x) 



dx itkc^S (x — X1X2) A(x\) B (X2) 



(9) 
(10) 

(11) 



The polarized coefficient functions and anomalous dimensions were calculated in [1-3] to next- 
to-leading order. 

The evolution equations may be rewritten changing the evolution variable log(Q 2 ) into 
a s using Eq. (§): 



d\og(Q 2 
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(12) 



We solve the evolution equations in MELLlN-iV space. For this purpose a Mellin- 
transformation 



M[f](N) = f 1 dxx N - 1 f(x), NeN 
Jo 



(13) 



of the equations is carried out under which the MELLIN convolution g) turns into the ordinary 
product. After the transformation was performed the argument A" is analytically continued to 
the complex plane. The principle way of solution of the Eqs. (0,f$D is described in the literature 
in detail, see 0, BTJ, which we summarize briefly. For the non-singlet case the solution to 
NLO is given by 
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(14) 



where a s = a s (Q 2 ),a = a s (Ql). The singlet solution reads 



Here the leading order evolution matrix L is given by 



L(a s ,a ,N) = e_(N) 



-r_ (AT) 



+ e. 



\a J 
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with eigenvalues of the LO singlet evolution matrix 

1 



tr(P (0) ) ± \tr(P {0) ) 2 - det 2 (P {0) ) 



(15) 



(16) 



(17) 



and the eigenvectors 
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The matrix Ui(N) is given by 

t/i(iV) = -e Pie - e + Pie_ 

with 
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(20) 



The input distributions Aq^ s (N, a ), A£(AT, a ) and AG(N, a ), see section 5, are evolved to 
the scale Q 2 , respectively to the coupling a s (Q 2 ). The inverse MELLlN-transform to x-space 
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is performed by a contour integral in the complex plane around all singularities, which can be 
written as 



1 



TT JO 



dz Im 



exp(i(j))x c{z) fcu[c{z 



(21) 



applying symmetry properties of the integrand. In practice an integral along the path c(z) = 
Ci+p [cos ((/>)+z sin(0)], with ci = 1.1, p > and (ft = (3/4)7r can be used. The upper bound on p is 
to be chosen by sufficient numerical convergence of the integral (|2l|) . The theoretical prediction 
/th^) for the respective observables depends on the parameters of the parton distributions 
chosen at the starting scale Ql and on Aq C d- These parameters are determined by a fit to the 
data using the x 2 ~ m ethod, see section 6. 

3 Scheme Invariant Evolution of Polarized Structure 
Functions 



In the foregoing section we investigated the QCD evolution of the polarized structure function 
gi(x,Q 2 ) at the level of twist-2 using the conventional picture of the QCD improved parton 
model. This requires choices of the non-perturbative partonic input distributions at a starting 
scale Qq. One outcome of the analysis is that some of the distributions are still very difficult to 
determine by a fit to the polarized deep inelastic data. This particularly applies to the polarized 
gluon and sea-quark densities. 

The scaling violations of deep inelastic structure functions, i.e. their Q 2 -behavior, are proper- 
ties of observables, and do thus not depend on particular representations, as scheme-dependent 
decompositions, which occur applying the parton model. Instead referring to the latter one 
might wish to eliminate the scheme dependence as occurring due to mass factorization all to- 
gether and describe the evolution of polarized structure functions based on observables only. As 
a consequence, no parton distributions would even emerge in this description, and choices of the 
parameterization of less constraint input densities can be avoided.^ 

Up to NLO the evolution of the structure function g\(x, Q 2 ) is described by one non- 
singlet (|7|) and the coupled singlet evolution equations (||). Instead of Eq. (|7|) one may use 
the scheme-invariant equation based on the non-singlet structure function gi S (x,Q 2 ) 



dt 



K m (a,x)®g^(x,Q 2 ) 
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Here the evolution variable is chosen as 



t 



Pi 



NS, 



X 



fa pNS, 



X) - A)Cq,l(x) 




(22) 



(23) 



(24) 



3 After having carried out the analysis projections onto the various parton-densities, including the gluon 
density, in whatever factorization scheme, the MS scheme or other schemes P^, ^7| , are possible including the 
respective experimental errors. 

4 Note that the splitting functions and coefficient functions are expressed in such a way that all evolution 
equations refer to logarithmic variations in Q 2 and a(Q 2 ) is used as the coupling constant. 
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For the singlet-evolution two observables have to be selected to form the factorization-scheme 
invariant evolution equations. In case of the polarized structure functions a natural choice is to 
use the flavor-singlet part of g±(x, Q 2 ) and its partial derivation w.r.t. t. f\ The singlet evolution 
equation reads (taking the Mellin transform) 
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(25) 



with Fa 
order by 



gi(N,Q 2 ), Fb = dgi(N,Q 2 )/dt. The physical evolution kernels are given in leading 
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In next-to-leading order they read 
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Here 7^ 0,1 ^ denotes the polarized deep-inelastic anomalous dimensions 
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5 Evolution equations of this type have been discussed before in Ref. B. In the unpolarized case one may 
consider the pair F^ix, Q 2 ), Fl(x, Q 2 ) as well, which has been considered in detail in Ref. |23| . For similar 
approaches, partly derived only for the small- a; domain, see ^8|. Corresponding evolution equations for time-like 
virtualities, which describe fragmentation in a factorization-scheme invariant way, were given in pa]. 
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and C^qM * s ^ ne MELLlN-transform of the Wilson coefficients 



°2,fc 



f dzz^Cftiz) . (33) 



For this combination in next-to-leading order the evolution depends on two evolution kernels 
only. 

The kernels K NS and Kg 1 are factorization scheme-invariant quantities. Their analytic struc- 
ture in z-space is difficult to obtain due to the inverse MELLlN-convolutions of inverse coefficient 
functions being required. Already in the case of the inverse of the leading order splitting function 



P gq {z) a rather complicated expression is found [29]. Therefore we perform all computations in 



MELLlN-iV space, where the physical evolution kernels are polynomials out of anomalous di- 
mensions and coefficient functions. Their analytic continuation to complex values of N can be 
performed using the respective representations of harmonic sums at high numerical precision, 
cf. p0||. In the unpolarized case approximate representations were also obtained in [Fill. 



The advantage of studying factorization-scheme invariant evolution equations both in the 
non-singlet and the singlet case is that the input distributions are observables. Although the 
present world statistics is too low, future high statistics measurements may provide accurate 
input densities gi(x, Qq) and dgi(x,Ql)/dt. In this case the only parameter to be measured 
analyzing the scaling violations of gi(x, Q 2 ) is Aqcd- At present the solution of the physical 
evolution equations in the polarized case cannot yet take full advantage of the method, since the 
respective observables are not yet measured well enough at typical input scales Qq. This different 
formulation, however, leads to an alternative view on the data in extracting Aqcd, as different 
input densities, gi(x, Qq) and dgi(x, Qfy/dt, are fitted as compared to AS(x, Qq) and AG(x, Qq) 
in the standard analysis. A comparison of the values of Aqcd obtained in both analyzes may 
indicate the stability of the determination of the QCD-parameter. 



4 Data 

The remarkable growth of experimental data on inclusive polarized deep inelastic scattering of 
leptons off nucleons over the last years allows to perform refined QCD analyzes of polarized 
structure functions in order to reveal the spin-dependent partonic structure of the nucleon. 
For the QCD analysis presented in the present paper the following data sets are used: the EMC 
proton data ||, the E142 neutron data ||, the HERMES neutron data ||, the E154 neutron data 
0, the SMC proton and deuteron data ||, the E143 proton and deuteron data |J, the HERMES 
proton data [10], the E155 deuteron data [11], and the E155 proton data [12 []. The number of 



the published data points above Q 2 = 1.0 GeV 2 for the different data sets are summarized in 
Table 1 for both the asymmetry data, i.e. g±/Fi or Ai, and data on gi together with the x— 
and Q 2 -ranges for the different experiments.^ There are 435 data points for asymmetry data, a 
number which exceeds the number of data points for g\ data by a factor of two. We therefore 
are utilizing the asymmetry data which are expected to give a better statistical accuracy. 

The QCD fits are performed on the polarized structure function gi(x,Q 2 ) which has to be 
evaluated from the asymmetry data. Experimentally cross section asymmetries for longitudinally 



6 Earlier data from Ref. [[13| are not considered. 

7 All corrections to the data are assumed to be carried out, including the QED radiative corrections 
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polarized lepton scattering off longitudinally polarized nucleons are measured 



A 



11 = ^TT^Tl • ( 34 ) 

The arrow-combination TT(TI) denotes parallel(anti-parallel) spin orientation of the lepton and 
the nucleon. The ratio of structure functions gi/F% and the longitudinal virtual photon asym- 
metry Ai are related to Au by 
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and the relation between gi/F\ and A\ is 
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Here A2 is the transverse virtual photon asymmetry and 

l-(l-y)e 
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l + eR(x,Q 2 ) 
4(l-y)- 



7V 



(35) 
(36) 

(37) 

(38) 
(39) 
(40) 

, 1 , ^ (41) 

e (! -y) 

D denotes the virtual photon depolarization factor, e, 7 and r\ are kinematic factors, y = Q 2 / (Sx) 
is the Bjorken variable, with S = M 2 + 2ME e , M the nucleon mass, E e the electron energy in the 
target rest frame, and R(x,Q 2 ) = (Jl/o't- F° r R(x,Q 2 ) the SLAC parameterization -R1990 is 
used by most of the experiments. At the time of the EMC experiment this parameterization was 
not available yet and R was assumed to be Q 2 independent. The SMC collaboration adopted a 
combination of -R1990 (for x > 0.12) and a parameterization derived by the NMC collaboration 
55H (for x < 0.12). In the E155 experiment a recent SLAC parameterization for R [36| was used. 



2y 2 + 4(1 
2Mx 
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7 



n 



ejy 



1 



The changes in the data caused by using the different R parameterizations are not significant 
and stay within the experimental errors^. 

and E143 [Pi and found 



The magnitude of A 2 has been measured by SMC |37||, E154 |38 



to be small. To a good approximation its contribution to gi/F\ and A\ which is being further 
suppressed by the kinematic factors 7 and rj can be neglected. Nevertheless, E143 and E154 have 
exploited their measurements, E155 has approximated this contribution by g^ w {x, Q 2 ) through 
the Wandzura-Wilczek expression p9| ,P| which is calculated from the measured structure function 
gi(x,Q 2 ) in the approximation that twist-2 contributions are dominant: 



S 2 WW (*,Q 2 



-gi(x,Q 2 ) + 



1 dz 



gi(z,Q 2 



(42) 



8 The EMC proton data, where the biggest impact is expected, change by a few percent only, see Ref. [[L5| . 

9 Note that this relation holds also in the presence of quark- |^o| and target-mass corrections [tlj . Related 
integral relations for twist-3 contributions and structure functions with electro-weak couplings were derived in 
Refs. Eol E2|. Recently these relations have been found to hold also for diffractive scattering, p3fl. 
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HERMES has accounted for the A 2 contribution to the proton data by using a fit based on 
existing measurements |TTJ. To obtain g±(x, Q 2 ) one has to multiply the ratio gi(x, Q 2 )/Fi(x, Q 2 ) 
with the unpolarized structure function Fi(x,Q 2 ) which can be calculated from the usually 
measured unpolarized structure function F 2 (x,Q 2 ) by 



;i+7 2 



2x(l + R(x,Q 2 



-F 2 (x,Q 2 



(43) 



For all data sets we are using the SLAC -Ri 99 o(a;, Q z 



and the NMC F 2 (x, Q 2 ) parameteriza- 
tion [H| to perform this calculation. 

The data sets used contain both statistical and systematic errors except the SMC data set 
which is given with statistical errors only. It is known that the systematic errors are partly 
correlated which would lead to an overestimation of the errors when added in quadrature with 
the statistical ones and hence to a reduction of the x 2 value in the minimization procedure. To 
treat all data sets on the same footing we decided to use the statistical errors only. However, 
we allow for a relative normalization shift between the different data sets within the normal- 
ization uncertainties quoted by the experiments. Thereby we are taking into account the main 
systematic uncertainties coming from the measurements of the luminosity and the beam and 
target polarization. The normalization shift for each data set enters as an additional term in the 
X 2 -expression which then reads 



X 
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(AN- 



+ E 



{Kg? 
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theor\2 

9i,j ) 



(Agf- 



data \ 2 



(44) 



where the sums run over all data sets and in each data set over all data points. The minimization 
of the x 2 value above to determine the best parameterization of the polarized parton distributions 
is done using the program MINUIT |45j . Only fits giving a positive definite covariance matrix at 
the end have been accepted in order to be able to calculate the fully correlated la error bands. 



5 Parameterizations of the Polarized Parton Distribu- 
tions and their Errors 

5.1 Parameterization of the Parton Densities 

The shape chosen for the parameterization of the polarized parton distributions in x-space at 
the input scale of Ql = 4.0 GeV 2 is 

xAft(x, Ql) = ViAx ai {l - x) lH (l + liX + Pi x*) . (45) 

The term x ai controls the low-x behavior of the parton densities and (1 — x) bi that at large 
values of x. The remainder polynomial factor accounts for the additional medium-x degrees of 
freedom. The normalization constants Ai 

Ar 1 = (l + 7i fl + B ^ b i + 1 )+ Pi B ( a * + \> b i + X ) ( 46 ) 

are chosen such that the rji are the first moments of Aq^x, Qq), rji = Jq dxAq^x, Ql). Here 
B(a, b) is the Euler Beta-function being related to the T-function by 
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L (a + o) 

In the present approach the QCD-evolution equations are solved in MELLlN-iV space as de- 
scribed in section |2|. The MELLlN-transform of the parton densities is performed and Mellin-N 
moments are calculated for complex arguments N : 



M[Af(x,Ql)](N) = f x N - l M{x,Ql)dx 

Jo 



Vi A i ( l + U TT, 1 + a ! ) B{N-1 + a h bt + 1) 
V N + di + hi J 

+(HB(N + a i -~ fc^ + l) ■ (4* 



As seen from Eq. (|45 ) there are five parameters for each parton distribution. To meet both the 
quality of the present data and the reliability of the fitting program the number of parameters 
has to be reduced. Assuming SU(3) flavor symmetry and a flavor symmetric sea one only has to 
derive one general polarized sea-quark distribution. The first moments of the polarized valence 
distributions Au v and Ad v can be fixed by the SU(3) parameters F and D as measured in 
neutron and hyperon /3-decays according to the relations : 

Vu v ~Vd v = F + D, (49) 
Vu v +Vd v = 3F-D. (50) 



A re-evaluation of F and D was performed in Ref. |15| on the basis of updated /3-decay constants 



1 49| leading to 



r) Uv = 0.926 ±0.014, (51) 
r) dv = -0.341 ± 0.018 . (52) 

Given the present accuracy of the data a number of parameters is set to zero, namely p Uv = 
Pd v = 0, 7 ? = pg = 0, and 7g = Pg = 0. This choice reduces the number of parameters to 
be fitted for each polarized parton density to three. In addition the QCD parameter Aq C d is 
to be determined. We allow for a relative normalization shift between the different data sets 
within the normalization uncertainties quoted by the experiments. These normalization shifts 
were fitted once and then fixed afterwards. 

The polarized parton densities to be determined are chosen to be 

Au v (x) = Au(x) — Au(x) , 
Ad v (x) = Ad(x) - Ad(x) , 

-AQ(x) = Aq(x) = Au(x) = Ad(x) = As(x) = As(x) , 
6 

and AG(x) . (53) 

Twelve parameters representing the parton densities are determined in the fit. These are : for 
Au v : a Uv , b Uv , and j Uv , for Ad v : a dv , b dv , and 7^, for AQ: rjq, a q , and b q , and for AG: i] G , 
da, and be- Note that rjg represents the first moment of the total quark sea. Starting off with 
these parameters the analysis shows, however, that the four parameters j Uv , 7^, bq, and be turn 
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out to have very large errors at Xmm- Altering them within the error range does not lead to 
a significant change of x 2 ; which shows that they are badly constrained by the inclusive data 
used in the analysis. We, therefore, fixed these parameters at their values obtained at the end 
of the fit and consider them as outer model-parameters. It is not expected that the small-x 
behavior of the polarized gluon and the sea-quarks is much different. To achieve this it turned 
out that the small-x slopes of the gluon and the sea-quarks are to be related like cic = aq + C, 
with C = 0.5 ... 1.0. We therefore decided to fix one of the parameter relative to the other. In 
fixing the high-x slopes be and bg we adopted a relation as derived from the unpolarized parton 
densities, namely bg/bcipol) = 6g/6c(unpol) (see e.g. Ref. PTfl ). Both relations together are 
suited to lead to positivity for AG and Aq. No positivity constraint was assumed for Au v and 
Ad v . After this the statistical measurement is only applied to the remaining seven parameters 
which describe the parton densities and to Aq C d- 



5.2 Error Calculation 



The evolved polarized parton densities and structure functions are linear functions of the input 



densities. Let f(x,Q 



2. 



7 \ k N 

L Ai=li 



be the evolved density at Q 2 depending on the parameters a 



Then its correlated error as given by Gaussian error propagation is 



Af{x,Q 2 



g(gw,) +i £( 



dcii da,j 



(Xj 



k 

l\i=l- 



(54) 



where C(ai, o 3 -) are the elements of the covariance matrix determined in the QCD analysis at 
the scale Qfo. The gradients df/ddi at this scale can be calculated analytically, except for Aq CD . 
Their value at Q 2 is calculated by evolution. The general form of the derivative of the Mellin 
moment M[/(a)](iV) w.r.t. the parameter a is 



dM[f(a)]{N) 
da 



F(a)xM[f(a)](N) 



(55) 



for complex values of N. Here we give the relevant expressions for F(a) as the normalized 
gradients w.r.t. the parameters being varied in the final fit, which are finally two parameters per 
distribution. For the polarized parton distributions Au v and Ad v one obtains 



F{ai) = tp(N - 1 + ai) - ij(N + en + h) 

n(h + 1) 



+ 



(N + ai + h)(N + ai + bi + 7i (JV - 1 + a^) 



- VK) + + h + 1) - ( -w^r^-ri — s) ' ( 56 ) 

\[Oi + k + l)(a; + bi + 1 + %a { ) J 

F(bi) = <%- + l)-V(A r + a i + 6 i ) 

/ ~fi(N -l + ai ) \ 



(N + ai + bi)(N + a l + b i + li (N -1 + a;)) 



^(fr. + l) + ^( a . + fr. + l) + ( -Tu^XI ■ i ■ S) ' ( 5? ) 

V {ai + bi + 1) [ai + bi + 1 + 7^) J 



10 For a discussion of positivity of parton densities see |5C 
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whereas for Aq and AG 



F((n) = ip(N -1 + di) -ip(N + ai + bi) -ilj(ai) +^(ai + bi + l) . 



(58) 
(59) 



Here ip(z) = dj dz (log T(z)) denotes the Euler -0-function. The gradients evolved in Mellin 
space are then transformed back to x space and are used according to Eq. (|54"1). To obtain 
the gradients for the error calculation of the polarized structure functions g\ and g™, which 
are composed of the polarized parton densities the expressions above, Eqs. (^- |59|) have to be 
multiplied with the corresponding Wilson coefficient functions. 

This yields the errors as far as the QCD parameter Aqcd is fixed and regarded as uncorrelated. 
The inclusion of the QCD parameter is performed easiest by numerical methods due to non- 
linear and iterative aspects in the calculation of a s (Q 2 , Aq CD ). The respective gradients are 
calculated performing the evolution both for A ± 8, with 5 < A 



df(x, Q 2 , A) f(x, Q 2 , A + 8)- f(x, Q 2 , A - 6) 



(60) 



OA 25 
using values for 5 ~ 10 MeV in the present analysis. 

Finally we also present the gradients at the input scale Qq in x-space for completeness, also 
w.r.t. the parameters being fixed in the final analysis : 



dAq t 1 dA qi 
= — Aq, , — — 

r]i ddi 



dA qi 
dh 



dr/i 
loefl — x) 



iog(x) - ?s 



1 dT s 

Tdbi. 



A qi 



dA Ql 



x 



with 



T 

dT 

da; 



dT 
db~. 



dT 

dT 

dpi 



dA qi 
dpi 



B(a i ,b i + 1)1 



X 2 



1 + 7jX + piX 2 T dp, 



yl+'JiX + PiXl 



}_dT\ 

TdJi) 



A qi 



x 



1 + CLi + h 

{en) - V>(oi + bi + 1)] B(a h bi - ■ \ ) ( I 



) +-fiB (ai + -,bi + l 



1 + + bi 



+ B{ ai ,bi + 1) 



+ 



1 3 



v 1 + at + 6i) 2 

[lf)(bi + 1) — i>(a,i + b; + I)] Z?(r, ; ,/,, + l)(l- 



p. fJ B ( a< + -,6j + 1 



(61) 
(62) 

(63) 

(64) 
(65) 
(66) 



1 + di + 6, 



V (l + Ai + fei) 2 
S(0i,6< + 1) 



^(6i + 1) - V(ai + ^ + -) 



S(ai,6i + 1) 



1 



Pifi ( aj + -A + 1 



1 + a,- + 



B[ai + -,bi + l 



(67) 
(68) 
(69) 



Both approaches have been used at the input scale and delivered the same error contours. 
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6 Results of the QCD Analysis 



6.1 Standard QCD Analysis 

In the fitting procedure we started with the 13 parameters selected, i.e. three parameters for each 
of the four polarized parton distribution and A QCD to be determined. For this set of parameters 
the sea-quark distribution was assumed to be described according to SU(3) flavor symmetry. 
This assumption is justified for an inclusive data analysis, which we aim on in this paper, given 
the present accuracy of the data, which may be refined in the forthcoming by including also semi- 
inclusive data, which become more and more precise. Analyzing the constraints of the different 
parameters at Xmin; ^ turns out that the deep-inelastic scattering data do not constrain the four 
parameters j Uv , 7^, bq, and bo sufficiently well, since rather large errors for these quantities are 
obtained. In the further procedure we fixed these parameters at their values obtained in the first 
minimization and chose the ratio bg/bc as in the unpolarized case [|6|. Let us note that the latter 



choice is one possible option. The lack of constraining power of the present data on the polarized 
parton densities has to be stressed, however. Only more precise data in this range can improve 
the situation in the future. In the fit we found some tendency towards a harder distribution for 
the gluon, although being not significant given the errors obtained for the large-x parameters. 
It also turns out that the small-x parameter for the gluon and the sea-quark distribution take 
essentially values which are shifted relative to each other as «g = as + c with c ~ 0.5 . . . 1, 
with some preference according to the value of xLm- This lead to values of c = 0.9, 0.6 in LO 
(ISET=1,2) and c = 0.9,0.5 for the NLO parameterizations (ISET = 3,4), cf. Appendix 1. 

The final minimization was carried out under the above conditions and determined the re- 
maining eight parameters with their la errors and the corresponding covariance matrix elements. 
Only fits ending with a positive-definite covariance matrix were accepted. The values and errors 
of these parameters along with those parameters fixed in the parameterization, cf. section 5, 
are summarized for both values of c in LO and NLO in table 2 for the parameters of the non- 
perturbative input distributions. The results on A QCD are discussed separately in section 6.2. 
The starting scale of the evolution was chosen as Ql = 4 GeV 2 . The covariance matrix elements 
for the LO and NLO fits, from which the parameterizations ISET=1. . .4 are derived, are given 
in tables 3-6. 

In Figures 1-4 the fitted parton distribution functions in leading and next-to-leading order 
for all sets of parameterizations and their la errors are presented at the starting scale Qq. The 
positivity bounds as e.g. obtained referring to the unpolarized distributions [H6] hold either for 



the central value itself or well within the present error bands. The current data constrain the up- 
valence distribution at best, followed by the down- valence, sea-quark, and gluon distributions. 



Our leading order results deviate from previous parameterizations |T^, [L6| somewhat at lower 
values of x. Except for the gluon distribution this deviation is much less in next-to-leading 
order, taking into account the experimental la errors. In the range x e [10~ 3 , 1] the polarized 
up-valence and gluon distributions are positive within errors, while the sea-quark and the down- 
valence densities are negative at the reference scale Ql = 4 GeV 2 , the latter except in a small 
range at very large x. 

In Figure 5 the NLO polarized gluon densities of ISET =3,4, are compared to the unpolarized 
gluon distribution of Ref. [fl6"|, to which the error of the gluon density as determined by the HI 



experiment [[IT] is overlaid symmetrically in the range x > 0.1. Both error contours illustrate 
the current margin of the positivity constraint for the gluon density, which is well covered. 
The polarized structure function xg\(x.Q 2 ) measured in the interval 3.0 GeV 2 < Q 2 < 
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5.0 GeV 2 , Fi gure 6, using the world asymmetry data is well described by our QCD NLO curve 
and the la error band. We also compare to corresponding representations of the parameteriza- 



tions [Uj], which are compatible within the present errors. 

Looking at the Q 2 dependence of the structure function gi(x,Q 2 ) in intervals of x gives 
insight to the scaling violations in the spin sector. As in the unpolarized case the presence of 
scaling violations are expected to manifest in a slope changing with x. The world proton data on 
gi(x, Q 2 ) have been plotted in such a way in Figure 7 and confronted with the QCD NLO curves 
of the present analysis and its la error bands. Corresponding curves of the parameterizations 



lq , |16f are also shown. Slight but non-significant differences between the different analyzes are 
observed in the intervals at low values of x. However, the data are well covered within the errors 
by all three analyzes. The current statistics in the low-x range is still rather low. 

In Figures 8-11 the scaling violations of the individual polarized momentum densities are 
depicted in the range x e [1(T 3 , 1], Q 2 e [1, 10 4 ] GeV 2 choosing the NLO distributions of ISET=3 
as an example. The up-valence distribution xAu v , Figure 8, evolves towards smaller values of 
x and the peak around x ~ 0.25 becomes more flat in the evolution from Q 2 = 1 GeV 2 to 
Q 2 = 10 4 GeV 2 . The distribution is positive within the errors. Statistically this distribution 
is constrained best among all others. The down- valence distribution xAd v , Figure 9, remains 
negative in the same range, although it is less constraint by the present data than the up- 
valence density. Also here the evolution is towards smaller values of x and structures at larger 
x flatten out. The momentum density of the polarized gluon xAG, Figure 10, is positive in 
the depicted range for Q 2 > 4 GeV 2 , but becomes slightly negative for smaller values of x at 
Q 2 ~ 1 GeV 2 . Also in this case the evolution moves the shape towards lower values of x and 
flattens the distribution. The error band becomes more uniformly. The sea-quark distribution 
xAq, Figure 11, is negative in the kinematic range shown for Q 2 < 10 GeV and remains negative 
within errors for x < 5 • 10-2 up to Q 2 = 10 4 GeV 2 , but chan ges sign for larger values of 
x. The minimum of the distribution at Q 2 = 1 GeV 2 around x ~ 0.1 moves to x ~ 0.01 at 
Q 2 = 10 4 GeV 2 . At the same time a maximum at x ~ 0.1 occurs. 

In the present analysis the structure functions gi ,n {x, Q 2 ) were parameterized in the twist-2 
approximation at next-to-leading order. As the data may contain higher twist terms as well it 
is of interest to look for the potential effect of these contributions. At present the size of twist-4 
contributions in deep-inelastic scattering data is widely unknown both in the unpolarized and 
polarized casePj Due to this we use the following phenomenological ansatze for a higher twist 
term : 



h(x,Q 2 ) = 1 + -4^(1 - , ("0) 

Q 
l 



h u (x,Q 2 ) = 1 + — A + Bx + Cx 2 , (71) 
Q 2 1 J 

which are used multiplicatively with gi(x, Q 2 ). We compare these fit results with the NLO 
parameterization ISET = 4 in Figure 12. Both ansatze yield results of similar size and deviate 
from the NLO curves at small values of Q 2 , however, they are fully consistent with the NLO 
results within the la error (hatched area). This shows that the present data do not contain 
significant higher twist contributions in the range Q 2 > 1 GeV and a NLO analysis can be 
carried out. This is also reflected by the observed logarithmic scaling violations, Figures 7, 12, 
and the measurement of the QCD scale Aqcd, to which we turn now. 

11 For a recent analysis of gluon induced contributions in the unpolarized case see Ref. fl54|. 
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6.2 A Q cd and cx s {M z ) 



In the QCD analysis we parameterized the strong coupling constant a s in terms of four massless 
flavors determining Aqq^ S . The NLO result fitting the asymmetry data, gi/F\ or A±, is 



A (4)MS 
ii QCD 

A (4)MS 
iv QCD 



235 ± 53 (stat) MeV, ISET = 3, 
240 ± 60 (stat) MeV, ISET = 4, 



(72) 
(73) 



identifying both the factorization and renormalization scales with Q 2 . The stability of the NLO 
result was investigated by changing both scales to values different of Q 2 . Since the present range 
in Q 2 probed in polarized deep inelastic scattering is still rather low we vary Q 2 only by factors 
of 2 and keep the other scale at Q 2 . The following variations are obtained: 



AA 



(4)MS 
QCD 



AA (4)MS 



+61 —61 

-47 (fac) + 114 (ren) MeV ' ISET = 3 ' 
_ 4° (fac) "^3 ( ren ) MeV ' ISET = 4 ' 



These results can be expressed in terms of a s (M§): 

a s {M z ) = 0.113 _ Q 004 (stat) _ Q 0Q4 (fac) _ Q 0Q5 (ren) 



ISET = 3, 



, , , , , 0.004 , , ^ +0.004 , +0.008 , , TOT _ A 
U. I I i (stat) _ Q 0Q4 (fac) _ Q Q06 (ren), ISET = 4 . 



<M 2 ) — _ 0>005 



Combining the errors these values 

a s (M: 



0.113 
0.114 



+0.010 
-0.008 , 

+0.010 



oc s (Ml _ _ oog 



(74) 
(75) 



can be compared with results from other QCD analyzes of polarized inclusive deep-inelastic 
scattering data 



E154 


[51]: 


<M 2 ) 


SMC 


@: 


a s (M 2 ) 


ABFR 


14 : 


a s (M 2 ) 



0.108-0.116, bad for > 0.120 , 
0.121 ± 0.002(stat) ± 0.006(syst + theor) 

120 +0 -° 04 (exp) +0 - 009 (theor) 
u - iAJ _ .005 lGXpj -0.006 ime ° rj ' 



and with the value of the current world average 

a s (M 2 z ) = 0.118 + 0.002 [53] . 



(76) 



(77) 



The results of the present analysis are consistent within la although our result and that of 
E154 |)lj] lead to a somewhat lower central value for a s (M z ). In the scheme-invariant analysis, 

which allows for an alternative view, a minor shift for AgS? down b Y 12 MeV 

was found leading; 
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to a central value of a s (M|) = 0.113 with statistical and the renormalization scale dependent 
errors of the same size as quoted above. Still more systematic effects have to be investigated 
in the future. Here we see a main point in studying further the unpolarized structure function 
Fi(x, Q 2 ), the denominator function in the expression for the spin asymmetry. In our analysis it 
was obtained from parameterizations of F2 [|4| and R |34j measurements which themselves are 
subject to systematic uncertainties. 



6.3 Scheme- Invariant Analysis 



A factorization-scheme invariant QCD analysis in next-to-leading order based on the observables 
gi(x, Q 2 ) and dgi(x, Q 2 )/dt for the proton has been performed, where t is the evolution variable 
defined in Eq. see section |3|. If compared to the standard analysis these two observables take 
the roles of AS(x, Q 2 ) and AG(x, Q 2 ). Such an analysis has, in principle, the advantage of direct 
experimental control over the input densities since they are measurable quantities. In this way no 
ansatz for AG is necessary and the only free parameter to be determined is Aqcd- Unfortunately, 
the quality of the present data does not yet allow an experimental determination of the slope 
dgi(x, Q 2 )/dt accurate enough to be used as input density. In the present analysis this slope has 
been derived fitting the world data as described in section 3. In the singlet analysis the initial 
distributions are gf(x,Ql) and dgf(x,Q 2 t )/dt. The latter quantity is depicted in Figure 13 as 
a function of x together with the corresponding slopes at higher values of Q 2 . This slope has a 
rather involved shape and it requires precise data in the range of the input scale Ql on gi(x, Q 2 ) 
to determine the slope experimentally in the future. The QCD analysis performed lead to a 
downward shift of 12 MeV for Aqcd which yields a similar result for a s (M§) as obtained in the 
standard analysis, see section 6.2. 



7 Moments of Polarized Parton Distributions 



In recent lattice simulations [23-25] low moments for the polarized parton densities Au v ,Ad v , 
and Au — Ad were determined. These moments are 



(x n+u 



>A/ = f x n+1 [Af{x) + (-l) (n+1) A/(x)] dx = M[Af{x)]{n) + (-l) (n+1) M[A/(x)](n) ,(78) 

J 



for n = —1,0, 1. The results of the simulation are presented at a scale fx 2 = l/a 2 ~ 4 GeV 2 , 
where a denotes the lattice spacing. The values of these lattice moments may be compared with 
the moments being obtained in the present analysis. Moreover we also present the respective 
moments for the sea-quark and gluon distributions and include their la error. Moments of 
structure functions or parton distributions based on a data analysis are only constrained in the 
range of the respective measurements. No definite statement can be made on the behavior of 
the (non-perturbative) parton distributions outside the range being explored]^] Reliably precise 



estimates neither for the range of very low [^0j nor large x are available currently. Due to this 



the value and error of the moments given below is that derived from the current experimental 
constraints. Nonetheless we study which contribution of the moments is determined by the range 
in which measurements exist. 



12 Note that before 1992 the strong growth of the structure function F%{x, Q 2 ), later found even at very low scales 
Q 2 ~ 2 GeV 2 , came as a surprise J55| . Likewise, at present the unpolarized gluon and sea-quark distributions at 
large values of x are widely unknown. 
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In Table 7 the moments of the NLO parton densities and some of their combinations, being 
derived in the present analysis, and their la errors are given. We also compare to the respective 
values obtained integrating only over the kinematic range in which data are measured and quote 
the value below and above the kinematic range for the moments of the NLO parton densities 
ISET=3,4. The determination of the error of the moments requires in general a correlated 
error propagation through the evolution equations, if the reference scale Q 2 is not the input 
scale Qq. As the lattice calculations provide measurements for a series of moments at the 
starting scale chosen in the present analysis the formulae given in section 5.2. apply directly. 
Whereas the estimated uncertainty for the distributions I SET = 3,4 for the unmeasured large-x 
range x > 0.85 is small, it is quite significant for the lowest moment for the small-x range, 
x < 0.02. For the up-quark distributions and their combinations the present correlated la 
errors are about 10% for the lowest moments, and reach about 30% for the down quark, and 
25% for the sea quarks. Depending on the parameterization the lowest moment of the polarized 
gluon distribution is found to have an error of 50% and more. Towards higher moments the 
errors grow, but the small-x uncertainties become less significant. The first three moments of 
the distributions Au v (x), Ad v (x) and Au(x) — Ad(x) can be compared with lattice results [23- 
25], with some caution. Only for a part of these values continuum extrapolations were yet 
performed. Still further study is required to safely determine and quantify systematic errors. 
While for the distribution Au — Ad the moment for n = — 1 comes out to be too small in the 
lattice measurements, it is somewhat too large for the moments n — 0, 1. We quote for the lowest 
moment an error, which is obtained due to all correlated parameters in the fit. On the other 
hand, the central value given is the axial charge g^, which is much more precisely determined by 
other measurements. It is yet difficult to see a systematic trend in the lattice values comparing 
to the moments found in the present analysis. Given the fact, that higher moments are more 
difficult to measure the agreement of the first two moments within the errors determined in the 
present analysis is quite good. 

We also compare the first moments in NLO with other recent analyzes, see table 8. Although 
some of the numbers, particularly for the gluon and sea-quark distributions look different they 
agree perfectly within the la errors derived in our analysis. 



8 Conclusions 

We have performed a QCD analysis of the inclusive polarized deep-inelastic charged lepton- 
nucleon scattering world data to next-to-leading order and derived parameterizations of polar- 
ized parton distributions at a starting scale Ql together with the QCD-scale Aq CD . The analysis 
was performed using the x 2 -method to determine the parameters of the problem in a fit to the 
data. A new aspect in comparison with previous analyzes is that we determine also the fully 
correlated errors of the parton densities and the QCD scale in leading and next-to-leading order. 
Due to the fact that not all shape parameters of the parton densities can be measured at suffi- 
cient accuracy using the present data, we derived two sets of parameterizations, which mainly 
differ in the low-x behavior of the gluon densities. Detailed comparisons were performed to the 
results obtained in other recent prameterizations [15", The previous results are widely com- 



patible with the present parameterizations within the current la error bands. Since we used the 
MELLlN-method to solve the evolution equations the Gaussian error propagation of the param- 
eters of the input densities through the evolution was possible in analytic form. Both the central 
values and the la errors of the parton densities are made available in form of a numerical pa- 
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rameterizations in the kinematic range x e [10~ 9 , 1], Q 2 e [1, 10 6 ] GeV 2 . These distributions can 
be used in the numerical calculations for polarized high-energy scattering processes at hadron- 
and ep-colliders. Moreover, due to the fact that parameterizations of the errors are available, 
error estimates of these quantities are possible w.r.t. the present knowledge of parton densities. 
These parameterizations are available as fast FORTRAN-routines which makes their application 
possible in Monte Carlo simulations. 

The current experimental data allow to measure the QCD scale Aq C d with a statistical 
error of SA = ±60MeV, and Sa s (M^) ± 0.004. Since at present only a next-to-leading order 
analysis can be carried out the variations in the renormalization and factorization scales induce 
yet large systematic errors. Combining all errors one obtains for the two scenarios (ISET = 3,4) 

a s (M|) = 0.113 ± q QQg and a s (M§) = 0.114 ± qqqq , respectively. The theoretical error due 

to the present renormalization and factorization scale uncertainties can be further reduced in 
a three-loop analysis, for which the anomalous dimensions have still to be calculated. We also 
performed for the first time a factorization scheme-invariant QCD analysis based on g\(x, Q 2 ) 
and dgi (x, Q 2 )/d log Q 2 which lead to a small shift of 12MeV in Aqcd only. This novel way 
of analysis which is based on the scaling violations of observables directly may show its full 
strengths in later analyzes based on the even higher statistics of future experiments. In this 
analysis factorization scale uncertainties do not occur. 

The results of the present analysis may be compared to recent lattice results calculating 
the first few moments of the distributions Au v (x, Q 2 ), Ad v (x,Q 2 ), and A(u v — d v )(x,Q 2 ) and 
their respective errors. We also present the respective moments of Aq(x, Q 2 ) and AG(x, Q 2 ) 
for which lattice results do not exist at present. Both on the side of the lattice measurements 
and the extraction of the distribution functions from deep inelastic scattering data the errors 
improved during recent years and the values became closer. However, more work has yet to 
be done in the future on systematic effects and even more precise experimental data would be 
welcome to perform an essential test of QCD also in this field at a higher precision than possible 
at present. 
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9 Appendix: The FORTRAN-code for the parton densities 
and their errors 



A fast FORTRAN program is available to represent the polarized parton densities xAu v (x, Q 2 ), 
xAd v (x, Q 2 ), xAG(x, Q 2 ), and xAq(x, Q 2 ) and the polarized structure functions xgi(x, Q 2 ) and 
xg™(x, Q 2 ) in leading and next-to-leading order in the MS-scheme together with the parameter- 
izations of their la errors. The following ranges in x and Q 2 are covered: 

1(T 9 < x < 1 , 1 GeV 2 < Q 2 < 10 6 GeV 2 . 

The polarized distributions are the result of a fit to the world data on spin asymmetries, i.e. 
A\' n ' d or gi/F[ ,n ' d , as described in the paper. The SUBROUTINE PPDF returns the values of the 
polarized distributions, always multiplied with x, at a given point in x and Q 2 by interpolating 
the data on specified grids. The interpolation in x is done by cubic splines and in Q 2 by a linear 
interpolation in log (Q 2 )- 

The parton distributions are evaluated by 

SUBROUTINE PPDF (ISET, X, Q2, UV, DUV, DV, DDV, GL, DGL, SEA, DSEA, 

G1P,DG1P,G1N,DG1N) 

All non-integer variables are of the type REAL*8. The calling routine has to contain the 
COMMON/INTINI/ I INI. Before the first call to SUBROUTINE PPDF the initialization is set by 
IINI = 0. The values of the parameter ISET are: 



ISET = 


1 


LO, 


a G 


= a s + 0.9 


ISET = 


2 


LO, 


a G 


= a s + 0.6 


ISET = 


3 


NLO, 




= a s + 1.0 


ISET = 


4 


NLO, 


a G 


= a s + 0.5 



The parameters X, Q 2 /GeV 2 are x and Q 2 . The momentum densities of the polarized up- and 
down valence quarks, gluons and the sea quarks are UV, DV, GL, SEA, with SEA = xAu s = 
xAd s = xAu = xAd = xAs = xAs. Correspondingly, DUV is the la error of UV etc. and G1P 
and GIN are the values of the electromagnetic structure functions g\ and g™. 
The program can be received on request via e-mail to Johannes.Bluemlein@desy.de or 



Helmut.Boettcher@desy.de or from [http : // www-zeuthen . desy . de/ ~hboett/ppdf . uu . gz 



13 We thank S. Kumano and M. Miyama of the AAC-collaboration for allowing us to use their interpolation 
routines. 
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Tables 



Experiment 


x-range 


Q 2 -range 
[GeV 2 ] 


number of data points 


Ref. 


gi/Fx or A 1 


9i 


E143(p) 
HERMES (p) 
E155(p) 
SMC(p) 
EMC(p) 


0.027 - 0.749 
0.028 - 0.660 
0.015 - 0.750 
0.005 - 0.480 
0.015 - 0.466 


1.17-9.52 
1.13 - 7.46 
1.22 - 34.72 
1.30 - 58.0 
3.50 - 29.5 


82 
39 
24 
59 
10 


28 
39 
24 
12 
10 


1 

I 


E 

1C 

1 
1 


1 

1 


proton 






214 


113 




E143(d) 
E155(d) 
SMC(d) 


0.027 - 0.749 
0.015 - 0.750 
0.005 - 0.479 


1.17-9.52 
1.22 - 34.79 
1.30 - 54.8 


82 
24 
65 


28 
24 
12 




i 

11 

i 




deuteron 






171 


64 




E142(n) 
HERMES (n) 
E154(n) 


0.035 - 0.466 
0.033 - 0.464 
0.017-0.564 


1.10 - 5.50 
1.22 - 5.25 
1.20 - 15.0 


30 
9 
11 


8 
9 

17 


§ 
I 

w\ 


51| 


neutron 






50 


34 




total 






435 


211 





Table 1: Published data points above Q 2 = 1.0 GeV 2 . 



20 





Scenario 1 


Scenario 2 




LO 


NLO 


LO 


NLO 




ISET=1 


ISET=3 


ISET=2 


ISET=4 




value 


error 


value 


error 


value 


error 


value 


error 


A$ CD , MeV 


203 


120 


235 


53 


195 


143 


240 


60 




0.926 


fixed 


0.926 


fixed 


0.926 


fixed 


0.926 


fixed 




0.197 


0.013 


0.294 


0.035 


0.199 


0.013 


0.271 


0.029 




2.403 


0.107 


3.167 


0.212 


2.416 


0.107 


3.070 


0.175 




21.34 


fixed 


27.22 


fixed 


21.34 


fixed 


27.22 


fixed 


Vd v 


-0.341 


fixed 


-0.341 


fixed 


-0.341 


fixed 


-0.341 


fixed 




0.190 


0.049 


0.254 


0.111 


0.182 


0.046 


0.325 


0.125 


bd v 


3.240 


0.884 


3.420 


1.332 


3.209 


0.895 


3.925 


1.129 


ld v (*) 


30.80 


fixed 


19.06 


fixed 


30.80 


fixed 


19.06 


fixed 


?]sea 


-0.353 


0.054 


-0.447 


0.082 


-0.314 


0.053 


-0.435 


0.061 


^sea 


0.367 


0.048 


0.424 


0.062 


0.428 


0.055 


0.285 


0.048 


b S ea(*) 


8.51 


fixed 


8.93 


fixed 


8.51 


fixed 


8.93 


fixed 


Vg 


1.281 


0.816 


1.026 


0.554 


1.043 


0.938 


0.931 


0.679 


ac 


a sea + 0.9 


a sea + 1.0 


a sea + 0.6 


a sea + 0.5 


M*) 


5.91 


fixed 


5.51 


fixed 


5.91 


fixed 


5.51 


fixed 


X 2 / NDF 


1.02 


0.90 


1.04 


0.93 



Table 2: Parameter values in LO and NLO (MS) of the 7 + 1 parameter fit based on the world 
asymmetry data for both scenarios. The (*) marks those parameters which were fixed after the first 
minimization since the present data do not constrain these parameters well enough (see text). 
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LO 




a (4) 
A QCD 




b Uv 




bd v 


Vsea 


O'sea 


Vg 


A (4) 


1.43E-2 
















a u v 


-2.05E-5 


1.80E-4 














b Uv 


-9.07E-5 


3.91E-4 


1.15E-2 












a d v 


1.10E-4 


1.03E-5 


-2.40E-3 


2.43E-3 










b d v 


-4.65E-5 


-7.92E-3 


-6.86E-3 


5.48E-3 


7.82E-01 








Vsea 


1.02E-4 


-4.46E-4 


-2.84E-3 


9.85E-4 


2.82E-2 


2.94E-3 






&sea 


-4.31E-5 


1.58E-4 


1.33E-3 


-5.96E-4 


-9.32E-3 


-2.58E-4 


2.29E-3 




Vg 


-1.03E-3 


2.02E-3 


1.58E-2 


-2.78E-3 


-1.61E-1 


-1.59E-2 


9.56E-3 


6.65E-1 



Table 3: The covariance matrix for the 7+1 parameter LO fit (scenario 1, ISET = 1) based on the 
world asymmetry data. 



NLO 




a(4) 
1V QCD 


a Uv 


b Uv 


ad v 


bd v 


Vsea 


&sea 


Vg 


A(4) 
^QCD 


2.81E-3 
















a Uv 


2.71E-5 


1.22E-3 














b Uv 


-1.30E-4 


5.10E-3 


4.50E-2 












ad v 


-3.35E-4 


-5.17E-4 


-3.23E-3 


1.23E-2 










bd v 


-6.22E-4 


-1.27E-2 


4.65E-2 


8.29E-2 


1.78E-0 








Vsea 


-5.30E-5 


-2.13E-3 


-1.12E-2 


5.19E-3 


4.74E-2 


6.77E-3 






^sea 


-4.85E-6 


9.07E-4 


4.49E-3 


-3.78E-3 


-2.98E-2 


-2.39E-3 


3.82E-3 




Vg 


4.03E-4 


1.41E-2 


6.71E-2 


-3.07E-2 


-2.22E-1 


-3.78E-2 


1.90E-2 


3.07E-1 



Table 4: The covariance matrix for the 7+1 parameter NLO fit (scenario 1, ISET = 3) based on 
the world asymmetry data. 
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LO 




a (4) 
A QCD 




b Uv 




bd v 


Vsea 


O'sea 


Vg 


A (4) 


2.05E-2 
















a u v 


1.25E-4 


1.76E-4 














b Uv 


-1.37E-3 


3.42E-4 


1.15E-2 












a d v 


-5.57E-4 


1.89E-6 


-2.19E-3 


2.14E-3 










b d v 


9.62E-3 


-7.84E-3 


-1.29E-3 


4.74E-3 


8.00E-01 








Vsea 


4.37E-4 


-4.34E-4 


-2.88E-3 


1.01E-3 


2.73E-2 


2.86E-3 






&sea 


-5.78E-4 


6.91E-5 


4.56E-4 


-3.73E-4 


-2.16E-3 


3.98E-4 


3.06E-3 




Vg 


9.43E-4 


1.91E-3 


1.54E-2 


-1.85E-3 


-1.55E-1 


-1.38E-2 


2.97E-3 


8.80E-1 



Table 5: The covariance matrix for the 7+1 parameter LO fit (scenario 2, ISET = 2) based on the 
world asymmetry data. 



NLO 




a(4) 
1V QCD 


a Uv 


b Uv 


ad v 


bd v 


Vsea 


&sea 


Vg 


A(4) 
^QCD 


3.55E-3 
















a Uv 


3.09E-4 


8.50E-4 














b Uv 


1.08E-3 


3.88E-3 


3.06E-2 












ad v 


1.88E-3 


2.66E-4 


-8.80E-4 


1.56E-2 










bd v 


1.74E-2 


-7.57E-3 


-1.34E-2 


9.03E-2 


1.27E-0 








Vsea 


5.77E-4 


-9.03E-4 


-5.32E-3 


3.26E-3 


3.84E-2 


3.74E-3 






^sea 


-5.33E-4 


6.81E-4 


3.68E-3 


-2.28E-3 


-2.56E-2 


-9.06E-4 


2.34E-3 




Vg 


-8.19E-3 


1.17E-2 


5.61E-2 


-3.95E-2 


-3.74E-1 


-2.37E-2 


1.51E-2 


4.61E-1 



Table 6: The covariance matrix for the 7+1 parameter NLO fit (scenario 2, ISET=4) based on the 
world asymmetry data. 
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I SET 


=3 


I SET 


=4 


lattice results 




n 


value 


value out 
of range 


value 


value out 
of range 


QCDSF 


LHPC/ 
SESAM 


Au v 


-1 


0.926 ± 0.071 


0.191|4E-4 


0.926 ±0.062 


0.207|4E-4 


0.889(29) 


0.860(69) 







0.163 ±0.014 


0.001|3E-4 


0.160 ±0.011 


0.001|4E-4 


0.198(8) 


0.242(22) 




1 


0.055 ± 0.006 


lE-5|3E-4 


0.055 ±0.005 


lE-5|3E-4 


0.041(9) 


0.116(42) 




2 


0.024 ± 0.003 


0|3E-4 


0.024 ± 0.003 


0|3E-4 






Ad v 


-1 


-0.341 ±0.123 


-0.104|-7E-5 


-0.341 ±0.103 


-0.086|-3E-5 


-0.236(27) 


-0.171(43) 







-0.047 ±0.021 


-5E-4|-6E-5 


-0.049 ±0.017 


-5E-4|-3E-5 


-0.048(3) 


-0.029(13) 




1 


-0.015 ±0.009 


-lE-5|-5E-5 


-0.015 ±0.007 


-lE-5|-2E-5 


-0.028(2) 


0.001(25) 




2 


-0.006 ± 0.005 


0|-5E-5 


-0.006 ± 0.004 


0|-2E-5 






Au-Ad 


-1 


1.267 ±0.142 


0.295|5E-4 


1.267 ±0.121 


0.293|5E-4 


1.14(3) 


1.031(81) 







0.210 ±0.025 


0.001|4E-4 


0.209 ±0.021 


0.001|4E-4 


0.246(9) 


0.271(25) 




1 


0.070 ±0.011 


2E-5|4E-4 


0.069 ± 0.009 


2E-5|4E-4 


0.069(9) 


0.115(49) 




2 


0.030 ± 0.006 


0|3E-4 


0.030 ± 0.004 


0|3E-4 






Au 


-1 


0.851 ± 0.075 


0.152|4E-4 


0.854 ± 0.066 


0.158|4E-4 











0.160 ±0.014 


8E-4|3E-4 


0.158 ±0.012 


8E-4|4E-4 








1 


0.055 ± 0.006 


lE-5|3E-4 


0.055 ±0.005 


lE-5|3E-4 








2 


0.024 ± 0.003 


0|3E-4 


0.024 ± 0.003 


0|3E-4 






Ad 


-1 


-0.415 ±0.124 


-0.144|-7E-5 


-0.413 ±0.104 


-0.135|-3E-5 











-0.050 ±0.022 


-7E-4|-6E-5 


-0.051 ±0.017 


-7E-4|-3E-5 








1 


-0.015 ±0.009 


-lE-5|-5E-5 


-0.015 ±0.007 


-lE-5|-2E-5 








2 


-0.006 ± 0.005 


0|-5E-5 


-0.006 ± 0.004 


0|-2E-5 






Aq 


-1 


-0.074 ±0.017 


-0.04|0 


-0.072 ±0.015 


-0.048|0 











-0.003 ± 0.001 


-2E-4|0 


-0.002 ± 5E-4 


-2E-4|0 








1 


-4E-4± 1E-4 


0|0 


-2E-4±6E-5 


0|0 








2 


-8E-5 ± 2E-5 


QjO 


-4E-5 ± 1E-5 


QjO 






AG 


-1 


1.026 ± 0.549 


0.04|lE-5 


0.931 ±0.669 


0.191|0 











0.184 ±0.103 


5E-4|lE-5 


0.100 ±0.075 


0.002|0 








1 


0.050 ± 0.028 


1E-5|1E~5 


0.022 ±0.017 


2E-5|0 








2 


0.017 ±0.010 


0|lE-5 


0.006 ± 0.005 


0|0 







Table 7: Moments of the NLO parton densities and their combinations for the parameterizations 
ISET=3,4 at Q 2 = 4 GeV 2 . The value of the respective moment integrating only outside the ar- 
range in which currently deep-inelastic scattering data are measured, 0.02 < x < 0.85, are given for 
comparison (lowerjupper part). The errors are the la correlated errors derived in the present analysis 
form the polarization asymmetry world data. The values of corresponding lattice measurements, 
cf. j[26|] , are shown for comparison. For n — 0, 1 for the values of the QCDSF collaboration no 
continuum extrapolation was performed. 
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Distribution 


ISET=3 


ISET=4 


ABFR [ 


14 


1 


GRSV | 


16 




AAC | 


15 




Au v 

Ad v 

Au 

Ad 

Aq 

AG 


0.926 ±0.071 
-0.341 ±0.123 

0.851 ±0.075 
-0.415 ±0.124 
-0.074 ±0.017 

1.026 ±0.549 


0.926 ±0.062 
-0.341 ±0.103 

0.854 ±0.066 
-0.413 ±0.104 
-0.072 ±0.015 

0.931 ±0.669 


7] u = 0.692 
r] d = -0.418 

1.262 


0.9206 
-0.3409 

0.8593 
-0.4043 
-0.0625 

0.6828 


0.9278 
-0.3416 

0.8399 
-0.4295 
-0.0879 

0.8076 



Table 8: Comparison of the first moments of the polarized parton densities in NLO in the MS scheme 
at Q 2 = 4GeV 2 for different sets of recent parton parameterizations. For the ABFR-analysis [|14| the 
values r) u> d are the first moments of An + Au and Ad + Ad, respectively, and As + As = —0.081. 
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Figure 1: LO polarized parton distributions at the input scale Q\ = 4.0 GeV 2 , ISET = 1, (solid 



line) compared to results obtained by GRSV (dashed-dotted line) |16| and AAC (dashed line) 
||15|| . The shaded areas represent the fully correlated la error bands calculated by Gaussian error 



propagation. The dark dotted lines indicate the positivity bounds choosing the distributions J27 
for reference. 
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Figure 2: LO polarized parton distributions at the input scale Ql = 4.0 GeV 2 , ISET=2, (solid 
line) compared to results obtained by GRSV (dashed-dotted line) |L6| and AAC (dashed line) 
The shaded areas represent the fully correlated lex error bands calculated by Gaussian 
error propagation. The dark dotted lines correspond to the positivity bounds according to the 
parameterization |27 ], 
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Figure 3: NLO polarized parton distributions at the input scale Ql = 4.0 GeV 2 , ISET=3, (solid 
line) compared to results obtained by GRSV (dashed-dotted line) JIB] and AAC (dashed line) 
fT5| . The shaded areas represent the fully correlated la error bands calculated by Gaussian 
error propagation. The dark dotted lines correspond to the positivity bounds choosing \T7\ for 
reference. 



28 




Figure 4: NLO polarized parton distributions at the input scale Ql = 4.0 GeV 2 , ISET=4, (solid 
line) compared to results obtained by GRSV (dashed-dotted line) |L6| and AAC (dashed line) 
The shaded areas represent the fully correlated lex error bands calculated by Gaussian 
error propagation. The dark dotted lines indicate the positivity bound if reference is taken to 
the distributions [f27fl . 
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Figure 5: NLO polarized momentum distribution for the gluon at the input scale Q$ = 5.0 GeV 2 
for ISET=3 (BB1, dashed line) and ISET =4 (BB2, dash-dotted line) with la error bands shown 
(shaded areas). The solid line corresponds to the unpolarized distribution [271. To the latter 



we added the experimental error of the unpolarized gluon distribution as determined in the HI 
experiment ||7| (hatched area), see also [48]. 
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Figure 6: The structure function xg\ measured in the interval 3.0 GeV 2 < Q 2 < 5.0 GeV 2 
as function of x. The structure function g±(x, Q 2 ) has been derived from the world asymmetry 
data. The error bars shown are the statistical errors only. The distribution is well described 
by our QCD NLO curve ISET=3 (solid line) at Q 2 = 4.0 GeV 2 and its fully correlated la error 
band calculated by Gaussian error propagation (shaded area). Also shown are the QCD NLO 
curves at the same value of Q 2 obtained by AAC (dotted line) [OJ and GRSV (dashed line) |TB| 
for comparison. 



31 



10 

: f Q 2 ) + 



i T SMC 

8 h ■■■^ - r - ■ * = 0-0063 (+7.5) £MC 

. j— x = 0.0141 (+6.2) 

AA ,.. r - M |-t- x = 0.0245 (+5.2) 



j-— y" x = 0.0346 (+4.5) 

£755 " ■ ■■■■hrjm x = 0.0490 (+4.0) 
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Figure 7: The polarized structure function g\ as function of Q 2 in intervals of x. The error 
bars shown are the statistical and systematic uncertainties added in quadrature. The data are 
well described by our QCD NLO curves (solid lines), ISET=3, and its fully correlated lex error 
bands calculated by Gaussian error propagation (shaded area). The values of C(x) are given in 
parentheses. Also shown are the QCD NLO curves obtained by AAC (dashed lines) |T5| and 
GRSV (dashed-dotted lines) |L6| for comparison. 
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Figure 8: The polarized parton distribution xAu v , ISET=3, evolved up to values of Q 2 = 
10, 000 GeV 2 (solid lines) compared to results obtained by GRSV (dashed-dotted lines) |L6| and 
AAC (dashed lines) |15[]. The shaded areas represent the fully correlated la error bands from 
our analysis calculated by Gaussian error propagation. 
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Figure 9: The polarized parton distribution xAd v , ISET = 3, evolved up to values of Q 2 = 
10, 000 GeV 2 (solid lines) compared to results obtained by GRSV (dashed-dotted lines) |16| and 
AAC (dashed lines) JT5fl . The shaded areas represent the fully correlated la error bands from 
our analysis calculated by Gaussian error propagation. 
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Figure 10: The polarized parton distribution xAG, ISET = 3, evolved up to values of Q 2 = 
10, 000 GeV 2 (solid lines) compared to results obtained by GRSV (dashed-dotted lines) Jl6| and 
AAC (dashed lines) fl5|| . The shaded areas represent the fully correlated la error bands from 
our analysis calculated by Gaussian error propagation. 
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Figure 11: The polarized parton distribution xAq, ISET =3, evolved up to values of Q 2 = 
10, 000 GeV 2 (solid lines) compared to results obtained by GRSV (dashed-dotted lines) |16f and 
AAC (dashed lines) fT5|| . The shaded areas represent the fully correlated la error bands from 
our analysis calculated by Gaussian error propagation. 
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Figure 12: Model fit to potential power corrections in gi(x,Q 2 ) as extracted from the world 
polarization asymmetry data in the present analysis (see text). Dashed line: model I, Eq. (f?"0|); 
dotted line: model II, Eq. ([71]). The full lines correspond to the parameterization (ISET=4) in 
the present analysis, to which the corresponding power correction model induces a perturbation. 
The shaded area corresponds to the la correlated error. 
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Figure 13: The evolution of dxgf(x,Q 2 )/dt, the slope of the polarized structure function 
xg\ (x, Q 2 ) of the proton (singlet contribution) with respect to t — — 2/j3 \n(a s (Q 2 ) / a s (Ql)) 
p3fl . The slope was determined from a fit to g\{x, Q 2 ), see text. 
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